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Abstract 


This paper suggests both a methodology to characterize ultracapacitors and to model their electrical behaviour. Current levels, frequency intervals, 
and voltage ranges are adapted to ultracapacitors testing. Experimental data results in the determination of the ultracapacitors performances in 
terms of energy and power densities, the quantification of the capacitance dependence on voltage, and the modelling of the dynamic behaviour of 
the device. Then, an electric model is proposed taking into account the ultracapacitors characteristics and their future use as peak power source for 
hybrid and electric vehicles. After, the parameters identification procedure is explained. Finally, the model validation, both in frequency and time 
domains, proves the validity of this methodology and the performances of the proposed model. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Among available storage devices, ultracapacitors exhibit 
interesting performances for a potential use in hybrid and elec- 
tric vehicles (HEV). The capacitance of large cells, which is 
about thousands of Farads, allows the storage of Wh energy with 
voltage limited to 2.7 V. Thanks to their very low serial resis- 
tance (<1 mQ), they can operate at high levels of instantaneous 
power. Also, they can work at very low temperatures (down to 
—40°C). Additionally, they can support a significant number 
of charge/discharge cycles. In typical applications as in hybrid 
and electric vehicles, they can be used as peak power source 
in combination with a primary source like battery, fuel-cell and 
generator [1,2]. 

An analysis of the power requirements in HEV applications 
and data from literature [3,4] show that the current profiles which 
are expected to be applied to ultracapacitors are composed of 
frequent charge and discharge pulses. This is so, because the 
ultracapacitors are storage devices with finite energy which need 
to be recharged periodically and because urban driving cycles 
are discontinuous speed profiles made of repetitive stop-and- 
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go micro-cycles with a typical period around 1-2 min. These 
charge/discharge pulses are characterized by high current lev- 
els, up to 600 A, and a width from about tens of milliseconds to 
tens of seconds. Following these specificities, we have focused 
on the dynamic behaviour of ultracapacitors at high current rates. 
This study is essential for efficient use of these components and 
to set-up thermal and reliability studies which are key issues in 
HEV. As ultracapacitors are new devices, their testing methods 
at high current rates are not well defined. Therefore, we present 
different methods dedicated to the characterization and the mod- 
elling of ultracapacitors at high current rates and the analysis of 
the corresponding results. Then, we propose an ultracapacitor 
electrical model adapted to a HEV use. Finally, the parameters 
identification and the validation of this model on typical current 
profiles are presented and discussed. 


2. Ultracapacitors characterization methods 


In this section, different methods are proposed for ultracapac- 
itors characterization and modelling [5-7]. The corresponding 
tests result in a lot of data and parameters that are very 
helpful to understand the ultracapacitors behaviour in different 
frequency intervals and different windows of observation-time. 
To do so, in this section, we present these tests, their practical 
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realization and the corresponding results for a Maxwell 2600 F 
cell BCAP0010. Measurements were conducted using a 
0-20 V/+800 A Digatron power tester and an IM6 impedance 
analyser with a 0—4 V/+40 A PP240 power potentiostat. 


2.1. Constant power tests 


Constant power tests are commonly used for the char- 
acterization of electric storage devices (batteries, flywheel, 
ultracapacitors, . . .) and generators (fuel-cells for instance). For 
ultracapacitors, the experimental procedure consists in discharg- 
ing the element with constant power between two voltage limits 
of Unom and Unom/2 in order to measure the specific energy 
that this element can supply as function of discharge power. 
As detailed in [8], the corresponding results can be graphically 
reported in the well-known Ragone diagram where the x-axis 
and the y-axis correspond, respectively, to the specific power 
(kW kg~!) and to the specific energy (Whkg~!). Nowadays, 
commercially available ultracapacitors have a typical specific 
power of 4kW kg™! and a typical specific energy of 6 Whkg™!. 
With such performances, they fill the gap between batteries and 
conventional capacitors. 

In practice, to obtain the points corresponding to the maxi- 
mum of specific power, the tester has to draw out a high current 
level under low voltage. However, as the time required to obtain 
these points is small, this maximum current value generally 
exceeds the peak current reached with repetitive profiles, as 
for instance, in power cycling tests. According to the specifi- 
cations given in [5], Fig. 1 shows the Ragone plot obtained for a 
2600 F/2.5 V cell and with a specific power that has been limited 
by the maximum current capability of the tester (800 A). Then, 
the experimental points have been interpolated and extrapolated 
according to Eq. (1), where E is the energy, P the power, R the 
series resistance and C is the rated capacitance [9]. 


s (prn +U? RP (1) 
= — n —; = 
2 Ue nom 


This diagram is particularly of interest in order to evaluate 
the suitability of a given storage device for an application by 
considering, for instance, the ratio between the specific energy 
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Fig. 1. Ragone plot of a 2600 F cell (obtained with maximum current of 800 A). 


and the specific power. In case of ultracapacitors, these time 
constants correspond to diagonal lines in Fig. 1 and show that 
about 6kW kg! are available during 0.2 s which is the typical 
duration for the ignition of the internal combustion engine in 
stop-start operation in HEV. Then, between 4 and 12 s, more than 
1kWkg~! can be supplied that is interesting for mild-hybrid 
vehicles. 

Finally, these constant power tests confirm that ultracapac- 
itors are convenient for applications that require peak power 
pulses with short durations of about few seconds as in micro- 
and mild-hybrids vehicles. 


2.2. Constant current tests 


Constant current tests represent a basic and a widely used 
characterization method that is useful to determine the rated 
capacitance and the equivalent series resistance (ESR). Then, 
these parameters can be determined both for charge and dis- 
charge and for different current levels and ambient temperatures. 
Therefore, the characterization procedure we have investigated 
on a 2600 F/2.5 V cell is shown in Fig. 2. 

This method is based on a succession of cycles, each made 
of a constant current charge up to Unom and a constant current 
discharge down to Unom/2, separated by a rest period during 
which the open-circuit voltage (OCV) is measured. 

Fig. 3 illustrates the evolution of the ultracapacitor volt- 
age measured at current switch-off after 200 A constant current 
charge and discharge with an ambient temperature of 20 °C. This 
method allows the determination of charge and discharge resis- 
tances ESR,, ESRg and charge and discharge capacitances Ce, 
Ca. Their values are given in Table 1 for different current levels. 

The results presented in Table 1 show that the ESR and the 
capacitance decrease with the increasing current level for both 
charge and discharge. This behaviour can be explained by con- 
sidering the porous nature of the electrodes which allows to 
reach a high specific area of capacitive interface. Therefore, this 
interface can be modelled by a volume distribution of inter- 
connected resistors and capacitors. The higher is the time of 
charge and discharge, the higher are the accessed resistance 
and capacitance. As the charge and discharge time decreases 
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Fig. 2. Experimental voltage response to 200 A charge/rest/discharge current 
cycle. 
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Fig. 3. Open-circuit voltage vs. time following current switch-off after charge 
(top) and discharge (bottom). 


with increasing current level when voltage limits are fixed, the 
measured capacitance and resistance decrease. In addition, the 
kinetic behaviour which governs the evolution of the OCV has 
been analysed and interpreted by considering different electro- 
chemical mechanisms like charge recovery and self-discharge 
[10]. 

Then, from these results, we can notice that there is a 
difference between the charge and discharge ESR values. Con- 
sequently, this shows a voltage dependency of the ESR as it was 
measured from Vi (1.25 V) for ESRa and from Vy (2.5 V) for 
ESR. 

The previously defined and calculated parameters can be 
associated to a simple electric model, suitable for rough simu- 


Table 1 
ESR and capacitance for different current levels 
Current (A) 50 100 150 200 
Vr. 
ESR, = JA (m2) 0.53 0.49 0.48 0.47 
c 
Vg 
ESRg = i (mQ) 0.46 0.44 0.43 0.42 
d 
Ces kT (F) 2825 2813 2803 2795 
© OCVy — OCVL 
la- Ta 
Ca = = (F 2824 2810 2801 2793 
1 = OCV — ocv, ® 


lations of the dynamic behaviour of ultracapacitors. In addition, 
the repeatability of this experimental method allows to compare 
different types of ultracapacitors or to estimate the ageing. 

In order to develop an enhanced electric model, we have to 
take into account the voltage dependency of the capacitance 
[11]. In fact, values of Ce and Ca reported in Table 1 are average 
capacitances calculated on the whole investigated voltage range. 
So, let us consider a voltage dependent capacitance. In this case, 
a differential capacitance Caif¢ is defined by: 


dq 
Caite(v) = — (2) 
dv 
After a constant current test, the measured capacitance is the 
mean value of the differential capacitance in the voltage range 
where the test is done. 
I(t¢ — ti) 1 hi 
= = - Caite dv (3) 
1 Jv; 


VE — Vj ve — v 


C 


Now, we suppose that the differential capacitance can be 
approximated with a polynomial function of voltage. 


Caitt = Co + Civ +++ + Cno” (4) 


For a constant current test, the integral of the differential 
charge leads to: 


fea f candv= fCotCivt t Cd (5) 


Clin 2 
Aq = CoV — Vi) + SW? — VP) +- 


C 
+4 A "a E ver) (6) 


According to the results of constant current tests, we can 
model the ultracapacitor with a series resistance and a voltage 
dependent capacitance. Therefore, from Eq. (6), the coeffi- 
cients C; can be determined by fitting the charge variation 
versus the ultracapacitor voltage. Also, the use of the relation 
C=I/(AV/At) on small voltage intervals allows the determina- 
tion of the differential capacitance versus voltage. Subsequently, 
Fig. 4 shows the variation of the stored charge in Ampere seconds 
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Fig. 4. Cumulated charge (A s) vs. voltage. 
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Fig. 5. Capacitance vs. voltage. 


versus the estimated ultracapacitor open-circuit voltage on the 
intervals of [1.25 V; 2.5 V]. We observe that a second order poly- 
nomial function can be used to fit the cumulated charge. Even if 
the difference with a linear fit is not evident, the residues, which 
correspond to the difference between measured and calculated 
points, prove that the second order polynomial function is the 
best fit. Consequently, we deduce the differential capacitance as 
a linear function of voltage. 

Fig. 5 compares the two methods. The results obtained using 
the relation C=I/(AV/At) validate the linear dependency of 
the capacitance with voltage. Also, they confirm a symmetric 
behaviour since there is no difference between charge and dis- 
charge. Also, Fig. 5 shows that the slope of the capacitance 
versus voltage curve is about 350F V~! which represents a 
more significant dependency than the one due to current level 
(Table 1). So, this result confirms the importance to take this 
non-linearity into account. 


2.3. Cyclic voltammetry 


Cyclic voltammetry is another test method commonly used 
by electrochemists in order to investigate the properties of 
electrodes. However, it is extended to ultracapacitors as the 
electrode/electrolyte interface governs the energy transfer mech- 
anism. This test is useful to validate the capacitive behaviour, to 
investigate the symmetry between charge and discharge process, 
and to determine the ultracapacitors potentials limits. It consists 
of linear voltage sweep between two limit values V; and V2. 
If the ultracapacitor is modelled with a constant capacitor C in 
series with an internal resistor Rs, the current is expressed by: 


mzc% = 7 
= =| = (=) 7) 


In this equation, there are two terms, the first is permanent 
I = Cdv/dt, the second is transient and disappears after a few 
time-constants. When the capacitance depends on voltage and 
assuming very low serial resistance, the current is proportional 
to the capacitance by 7 = Caițtdv/dt given that the voltage drop 
due to the resistor can be neglected. 
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Fig. 6. Voltamogram measured with a sweep rate of 10 mV s~! ona 2600 F/2.5 V 
cell. 


Fig. 6 illustrates the voltamogram which is obtained by 
recording the current versus voltage with a sweep rate of 
10mVs~!. This figure shows that the current is constant 
between 0 and 0.3 V, i.e. the differential capacitance is constant 
within this voltage range. Then, the current increases as well as 
the differential capacitance. In the [1.25 V; 2.5 V] range, the lin- 
ear fit of the differential capacitance which has been determined 
from constant current tests, seems to be right. Also, the evolu- 
tion shape of the current versus voltage seems to be symmetric 
between charge and discharge that corresponds to a reversible 
behaviour, i.e. the differential capacitance in charge is the same 
as in discharge. 


2.4. Impedance spectroscopy tests 


The existing impedance analyzers allow the investigation of a 
wide frequency range from the Hz to the MHz. But, generally 
the output current of these apparatus is very low and limited to 
few hundred of mA. In order to measure very low impedance 
(<mQ), voltage variations (in the range of few pV) have to be fil- 
tered from a response signal with a large dc bias voltage. These 
requirements are very difficult to satisfy and the lower is the 
measured impedance the lower is the accuracy. For this rea- 
son, current boosters able to apply high ac currents excitations 
are used in order to ensure an accurate impedance measure- 
ment. In this study, impedance measurements are made with 
an electrochemical workstation associated with a + 40 A power 
booster. 

Firstly, Fig. 7 depicts the module and phase of the ultracapac- 
itor impedance versus frequency and Fig. 8 shows the real and 
imaginary parts versus frequency. The measurements have been 
done for a 2600 F/2.5 V ultracapacitor in the [1 mHz; 60 kHz] 
frequency range. The voltage stimulus was characterized by 
a 10mV ac peak value and has been applied around 2 V dc 
bias potential. Consequently, three frequency intervals can be 
distinguished: 


e Very low frequencies (<10 mHz): in this range, the impedance 
magnitude is relatively high (>5 m2). Then, as the phase 
angle is close to —90°, the ultracapacitor impedance is mainly 
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Fig. 7. Module and phase of the ultracapacitor impedance. 


capacitive. The particularity of this region is an increase of the 
real part conversely with frequency which traduces a typical 
charge redistribution phenomenon [12]. 

e Medium frequencies [10 mHz; 1 kHz]: we observe that the 
ultracapacitors behaviour changes from capacitive to induc- 
tive. This range is characterized by very low impedance which 
reaches a minimum value of 0.27 mQ for a null phase angle. 
Thus, the ultracapacitors are very convenient to be used as 
peak power source when the current profiles have their fre- 
quency spectrum in this interval. 

e High frequencies (>5 kHz): the impedance magnitude is rela- 
tively high (>2 mQ), the phase angle is close to +80° and we 
notice an inductive behaviour as the impedance real part is 
increased with the frequency. 


Then, we have focused on the impedance in the [10 mHz; 
1 kHz] frequency range which is more adapted to their use in 
HEV applications since the required power pulse will have their 
frequency spectrum in this interval. Fig. 9 shows this measured 
impedance in the Nyquist plan. 

At low frequencies measurement points follow a vertical line 
which is the typical Nyquist plot of a pure capacitance in series 
with a resistance named ESRgc. Then, for intermediate frequen- 
cies, the impedance real part decreases with frequency as well as 
the imaginary part. This behaviour corresponds to the diffusion 
in a porous electrode and can be approached by an equivalent 
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Fig. 8. Impedance real and imaginary parts vs. frequency. 
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Fig. 9. Nyquist diagram of the ultracapacitor impedance. 


RC transmission line. Finally, beyond a resonance frequency of 
about 200 Hz, the nature of the impedance changes and becomes 
inductive whereas the real part tends towards a minimal value, 
namely ESRprF. 

Several measurements have been done at various bias volt- 
ages in order to determine the impedance dependency on 
voltage. As illustrated in Fig. 10, the imaginary part of the 
impedance, in absolute value, increases when the dc bias poten- 
tial is reduced. This result traduces a capacitance increase as 
the voltage increases, as it was previously highlighted thanks 
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Fig. 10. Impedance imaginary part vs. frequency. 
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Fig. 11. Impedance real part vs. frequency. 


to constant current and cyclic voltammetry tests. Also, as illus- 
trated in Fig. 11, we can notice that the impedance real part at 
low frequency depends on voltage. 


3. Ultracapacitor electrical modelling 
3.1. Proposed electric model for ultracapacitors 


Different studies were devoted to ultracapacitor modelling 
{11,13-15]. An overview of these works combined with an anal- 
ysis of experimental data shows that the ultracapacitor electric 
model has to take into account many phenomena which are: 


e the ultracapacitor behaviour depends on many physical 
parameters basically voltage and temperature; 

e an equivalent RC transmission line behaviour that character- 
izes the ultracapacitors dynamic response especially in the 
[0.1 Hz; 10 Hz] frequency range. This behaviour is induced 
by the porous nature of the capacitive interface; 

e aresonance frequency (<200 Hz) that corresponds to a transi- 
tion between capacitive and inductive nature of the impedance 
and measured at the minimum of the impedance real part; 

e an increase of the impedance real part versus frequency and 
an inductive behaviour above the resonance frequency; 

e a charge redistribution phenomenon that occurs at low fre- 
quencies or for charge and discharge higher than | min. It was 
modelled by two RC branches characterized by long constant 
times compared with the constant time of the RC transmission 
line; 

e a self-discharge which can be modelled by a high parallel 
resistor called leakage resistor. 


Therefore, a general electrical model is proposed in Fig. 12. 
But, according to the specifications of the application we can 
neglect some phenomena in order to use the most suitable and 
the least complicated electric model. In this aim, the window 
of time-observation of the current and of the voltage that we 
need for the application is useful. Regarding HEV applications, 
as previously mentioned, we will focus on the frequency range 
from 10 mHz to 1 kHz. Therefore, the charge redistribution and 
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Fig. 12. General electric model of ultracapacitor. 


the self-discharge phenomena will be neglected. Also, the high 
frequency behaviour will not be modelled. So, the proposed elec- 
tric model is given in Fig. 13. It is made of a serial inductor 
L,, an access resistor Race which corresponds physically to the 
resistance of the device terminals, a capacitor Cacc and a non- 
linear transmission line with voltage dependent capacitors. The 
impedance real part dependency on voltage is neglected. Finally, 
the transmission line is approached by n RC branches. The higher 
is the order n, the better the dynamic behaviour is modelled. 


3.2. Model parameters extraction and validation 


The model parameters are identified in different steps using 
experimental data essentially those obtained from constant 
current tests and impedance spectroscopy. An identification pro- 
cedure programmed in Matlab leads to single combination of 
parameters. At first, we determine the value of the access resis- 
tance Racc which corresponds to the minimum of the impedance 
real part or the magnitude of the measured impedance at the reso- 
nance frequency. Then, for frequencies higher than the resonance 
one, we assume that the transmission line can be neglected. 
Therefore, the serial inductor Ls and the access capacitance can 
be extracted by fitting the imaginary part versus frequency. For 
the tested ultracapacitor, values of 200 F and 50 nH are obtained. 
Fig. 14 validates the used method by comparing measured points 
and model results. 

Finally, we determine the transmission line parameters R and 
C. The transmission line capacitance, C, is deduced thanks to 
the test of charge and discharge with constant currents. In this 
test, the measured capacitance is the sum of all the capacitance 
appearing in the electric model. If we neglect the phenomenon 
of charge redistribution, the measured capacitance will be the 
sum of Cacc and the capacity of the line of transmission C. For 


Fig. 13. Ultracapacitor electric model proposed for HEV applications. 
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20). 


The dependence of the capacitances with voltage is verified 


by simulating acyclic voltammetry. Fig. 18 compares the experi- 
mental and simulation results. A difference between simulations 
and experiments is noticed when the voltage is lower than 0.8 V; 
then the difference becomes very slight. Indeed, in the model 
where the ultracapacitor is charged and discharged between 1.25 


and 2.5 V. To enhance the electric model, 
ues should be obtained from cyclic voltammetry test. But in this 


and 20 (Fig. 17). Nevertheless, the user has the possibility to fix 
the capacitances values are obtained from constant current tests 


have been obtained for an order comprised between 4 (Fig. 16) 
the order arbitrary in the identification program. 


Fig. 17. Model validation in frequency domain: real part vs. frequency (trans- 
1 kHz]. Concerning the transmission line order, good results 


mission line order 


10 
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Fig. 14. Comparison of measured and simulated imaginary parts at high fre- 


quencies. 


study the model accuracy for voltage low than 1.25 V is not very 


interesting and so the obtained results are accepted. 


Fig. 15. Model validation in frequency domain: imaginary part vs. frequency. 


Finally, in order to validate the model with high current rate, 
we have compared the voltage response of the model to the 


the tested ultracapacitor, we obtain: 


experimental response obtained with a 2600 F/2.5 V cell with 
a specific current profile made of 200 A charge and discharge 


pulses. This comparison is illustrated in Fig. 19. By modifying 
the duration of the charge and discharge pulses, we can estimate 


C(v) [F] = 2100 + 360 x V = C = 1900 + 360 x V 


The value of the resistance R is determined by fitting the 
measured real part versus frequency to the impedance real part 
of the proposed electric circuit. For the tested device we obtain 


a value of R 


the precision of the model for different voltage variations. This 
result shows that, after 1000s, there is no global divergence 


0.54mQ. 


Figs. 15 and 16 validate the proposed model in the frequency 
domain and in the investigated frequency interval [10 mHz; 
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Fig. 18. Cyclic voltamogram: comparison between simulation and experimental 


results. 


Fig. 16. Model validation in frequency domain: real part vs. frequency (trans- 


mission line order = 4). 
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Fig. 19. Model validation with discontinuous current profile (200 A charge and 
discharge pulses). 


between measured and calculated points. Also, we can observe 
a good matching at the time scale of few seconds. 


4. Conclusion 


In this paper, ultracapacitors testing methods dedicated to 
characterization and modelling have been presented. Current 
levels, frequency intervals and voltage ranges have been defined 
according to the HEV requirements and the future use of ultraca- 
pacitors as peak power devices. Therefore, we have focused on 
the dynamic behaviour of ultracapacitors and we have noticed 
that they have a particular behaviour which is essentially capac- 
itive at low frequencies, equivalent to a transmission line at 
intermediate frequencies, and inductive at high frequencies. 
Also the capacitance dependence with voltage has been high- 
lighted. 

Then, taking into account the ultracapacitors characteristics 
and the specificities of current profiles, we have proposed a 


dynamic electric model. The model validation in both frequency 
and time domains has proved the interest of this methodology 
and the performances of the model. 

The investigation of the ultracapacitors electrical behaviour is 
a preliminary step before their integration in real applications. 
The model has to be completed with temperature dependent 
parameters. This study has to be completed with thermal and 
reliability studies as high current rates lead to self-heating and 
ageing. 
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